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Abstract 

For a second-order linear differential equation with two irregular singular points of rank 
three, multiple Laplace-type contour integral solutions are considered. An explicit formula 
in terms of the Stokes multipliers is derived for the characteristic exponent of the multi- 
plicative solutions. The Stokes multipliers are represented by converging series with terms 
for which limit formulas as well as more detailed asymptotic expansions are available. Here 
certain new, recursively known coefficients enter , which are closely related to but different 
from the coefficients of the formal solutions at one of the irregular singular points of the 
differential equation. The coefficients of the formal solutions then appear as finite sums 
over subsets of the new coefficients. As a by-product, the leading exponential terms of the 
asymptotic behaviour of the late coefficients of the formal solutions are given, and this is a 
concrete example of the structural results obtained by Immink in a more general setting. 
The formulas displayed in this paper are not of merely theoretical interest, but they also 
are complete in the sense that they could be (and have been) implemented for computing 
accurate numerical values of the characteristic exponent, although the computational load 
is not small and increases with the rank of the singular point under consideration. 

AMS classification: 34A20; 34A25; 34A30 
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1. Introduction 

Let us consider the differential equation 

6 6 

Z 2 f" + Zf'-[J2 D ^ Z ~ m + L2 +J2 B mZ m ]f{z) = 
rn=l m=l 

1 



(1.1) 



with the thirteen parameters Di, . . . , Dq, L, B±, . . . , Bq, which for simplicity of presenta- 
tion are assumed to be real. This differential equation has two irregular singular points, 
each of rank 3 (if Dq ^ and Bq ^ 0, respectively), at the origin and at infinity. Without 
loss of generality, one of the parameters except L could be set equal to 1 . We assume that 
Bq is positive, again for simplicity of presentation, and that L is not negative. 
At infinity, there are formal power series solutions 

oo 

= eMP(z))z- T{1) £ a n (l)z-\ (1.2) 

n=0 

oo 

fZl (z) = eM-P(z))z- Ti - 1] E *n(-l)z- n , (1.3) 

n=0 

where the various quantities are determined by 

P(z) = p 3 z 3 + p 2 z 2 + Pl z, (1.4) 
p 3 = \\/Bq~, P2 = j2B 5 /p 3 , pi = \{B A - 4p|)/p 3 , (1.5) 

1 Kps 6 KP3 

KG {1,-1} (1.7) 



and where the coefficients 



a„(«) = a n (1.8) 



are known recursively by 

a = 1 (1.9) 

and 

6Kp 3 na n = [-Akp 2 {t(k) + n - 2) +p\ - B 2 }a n -i 

+ [—2kpi(t(k) + n - |) - Si]a n _2 + (t(k) + n — 3 — L)(t(k) + n - 3 + L)a n _ 3 
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(1.10) 



^ ^ Dm&n—m — 'Si 



m=l 



n = 1, 2, . . . , (a_8 = a_7 = . . . = a_i = 0). 
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The formal solutions are asymptotic expansions as z — > oo in appropriate sectors of the 
complex plane. 

At the origin, there are analogical formal solutions, which may be obtained from those 
at infinity by the simultaneous replacements z «-> 1/z, B m <-> D m , m = 1, . . . , 6. This 
symmetry is our main reason for choosing just (1.1) as the standard form for a differential 
equation with two irregular singular points. 

In the ring-shaped region < \z\ < oo we have (convergent) Floquet solutions f^\z) 
and f(~ u '(z), where 

/<">(*) =Z" C n )z ^ ( L11 ) 

n= — oo 

which are linearly independent if 2u is not equal to an integer. Here the coefficients e4 
obey the recurrence relation 

6 6 

- B^n-m + (-L + uj + n)(L + u + n)<£ w > - £ D m c^ m = 0, (1.12) 

m=l m=l 

and we want to normalize them by choosing 

= 1. (1.13) 

The requirement that the power series converge determines, modulo 1, the possible values 
of the characteristic exponent (or circuit exponent or Floquet exponent) ui. The problem 
to compute the characteristic exponent appears also and is best known in the context of 
Hill's differential equation [22]. There are methods to compute the characteristic exponent 
numerically, which require, for ui = 0, the evaluation of the infinite determinant [20] 
associated with (1.12), or numerical integration of the differential equation along a suitable 
contour [21], or numerical solution of an eigen- value problem [14]. 

This paper developes an entirely different method for evaluating the characteristic ex- 
ponent. We obtain an explicit formula in terms of quantities which are essentially the 
Stokes multipliers, and these are given explicitly as convergent series, the terms of which 
are represented by asymptotic expansions. Here certain recursively known coefficients en- 
ter which are closely related to but different from the coefficients of the formal solutions 
(1.2)-(1.3). 

3 



Although some of the other authors concerned with irregular singular points of rank 
larger than one [1, 2, 4-6, 8-10, 12, 13, 15, 18, 19] consider the general case of arbitrary 
rank , we here prefer to restrict our attention to rank three. This is already general enough 
to give an impression of what can be expected in the case of even higher rank. On the other 
hand, it is still simple enough so that we can, for the relevant quantities, obtain explicit 
expressions which are not only theoretically interesting but can also be implemented (and 
have been implemented) for numerical evaluation. This work may be viewed as an attempt 
to extend [3] , which was useful for rank one, to the much more complicated case of higher 
rank. 

2. Laplace contour integral solutions 

We try to apply the classical method of multiple Laplace contour integral solutions [7] 
and write for a solution of (1.1) 



where a power factor with a still arbitrary parameter A has been included in view of later 
benefits. To derive the appropriate weight function v(s,ti,t 2 ) is a somewhat lengthy but 
not principally difficult procedure. We therefore give two lemmata stating the results of 
this procedure and postpone the proofs to a later section. 

Lemma 1. The weight function v(s, ti, t 2 ) has to be a solution of the partial differential 
equation 



9(s 2 - sl)(d A v/ds A ) + 12(s*i - s t 10 )(d 4 v/ds 3 dt 1 ) + 6(st 2 - s t 20 )(d 4 v/ds 3 dt 2 ) 
+ A{t\ - t 2 )(d 4 v/ds 2 dt 2 ) + {t\ - t 2 )(d 4 v/ds 2 dt 2 ) + 4(*i* 2 - t 10 t 20 )(d 4 v/ds 2 dt 1 dt 2 ) 
+ ([81 - 6X]s + [B 3 - 4t 10 t 20 ])(d 3 v/ds 3 ) + ([52 - 4A]*i + [B 2 - t 2 20 ])(d 3 v/ds 2 dt 1 ) 
+ ([25 - 2X]t 2 + B 1 )(d 3 v/ds 2 dt 2 ) + ([A - 12] 2 - L 2 )(d 2 v/ds 2 ) - D^v/dsdh) 
- D 2 (d 2 v/dsdt 2 ) + D 3 (dv/ds) + D^dv/dh) + D 5 (dv/dt 2 ) - D 6 v = 




(2.1) 



(2.2) 



with two finite singular points at (s,ti,t 2 ) = (ksq, Kt\o, Kt 2 o) where 



ho = Pi, 



ho = Pi, 



so =P3, 



(2.3) 
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and the contours for each variable have to satisfy the condition that a certain lengthy 
expression, bilinear in K and v or their partial derivatives, have the same value at the 
start and the end of the contour. 

Lemma 2. For k G {1,-1}, there are appropriate weight functions v = V(k; s, t\, t 2 ) 
which at the singular point (s,ti,t 2 ) = (ksq, utio, ^20) have the power series expansion 

V(K;s,ti,t 2 ) 

00 00 00 

= E E E M^rn,n 1 ,n 2 )(s- K sor {K)+m (t 1 - K t 1 o)-^- n Ht2-^ 20 )-^- n \ 

m=0 ni=0 ri2=0 

(2.4) 

where the A- coefficients satisfy a certain recurrence relation and the exponents are related 
by 

3//(k) - 2v x - vi = X + t(k) - 6. (2.5) 
An appropriate set of coefficients is 

A(k; m, ni, n 2 ) = T(-/j(k) - m)T{vi + rii)r(z/ 2 + n 2 )b(n; m, ni, n 2 ), (2.6) 

where the new coefficients b(K;m,ni,n 2 ) are given by the recurrence relation 
6ftSo(3m — 2ni — n 2 )b(K; m, n\,n 2 ) 
= [(3m — 2n\ — n 2 + t(k) — 3) 2 — L 2 ]b(n; m — 1, ni, n 2 ) 
+ [-AKt w (3m - 2m -n 2 + t(k) - 2) + t\ - B 2 ]b(n; m — 1, ni — 1, n 2 ) 

(2.7) 

+ [— 2Kt 20 (3m — 2ni — n 2 + t(k) — |) — I?i]6(k; m — 1, ni, n 2 — 1) 

— Dib(n; m — 2, ni — 1, 712) — D 2 b(K; m — 2,ni,n 2 — 1) — D 3 b(n; m — 2,n\, n 2 ) 

— D$b(K; m — 3, ni — 1, 712) — D^b^n; m — 2>,n\,n 2 — 1) — D^b(K; m — 3, ni, 712) 
witt t/ie initial conditions 

b{n; 0,0,0) = 1, 

m, ni, 77-2) = 2/ 3m — 2n± — n 2 = for (m,ni,n 2 ) ^ (0,0,0), (2.8) 
(6(k; m, ni, 71-2) = z/ m<0 or ni < or ri2 < 0) 

and u>irj/i 

6(/c;0,m,n 2 ) = /or (m, n 2 ) ^ (0, 0), (2.9) 
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b(K; m, tt-1, 712) = for n\ + n 2 > m (2-10) 

as a consequence. In addition, there are weight functions U(k; s,ti,t 2 ) which are analytic 
in s at the respective singular point, corresponding to h(k) = 0,1,2 (without any further 
relation such as (2.5)). 

In order to avoid unnecessary complications, we want to assume that the non-trivial 
exponent fi according to (2.5) is not equal to an integer. This can always be guarantied 
by a suitable choice of the still disposable parameter A. 

Since the exponents //(k), vi, v 2 are restricted only by (2.5) but otherwise arbitrary, 
there are other solutions of the partial differential equation (2.2 ) relevant as weight func- 
tions in our contour integrals (2.1). We may assume that v\ and v 2 are positive integers, 
preferentially 

vx = vi = 1, (2.11) 

but for the time being we want to keep v\ and v 2 in the formulas. If v\ and v 2 are increased 
by any positive integers q\ and q 2 , respectively, and simultaneously is decreased by 
(2/3)(/i + (1/3)^2, then (2.5) is still satisfied. We therefore have to consider, for 91,52 = 
0, 1,2,.. ., the set of solutions 

V(k; qi, q 2 ) = V(k; qi,q 2 ; s, t u t 2 ) 

00 m m 

= ^2 ^2 ^2 + qi + ?ll ) r ( z/ 2 + Q2 + n 2 )T(-/j(K) - §<?i - \q 2 - m) 

m=0 ni=0 n2=0 

x b(K;m,n 1 ,n 2 )(s - ^o) MW+ ^ 1+ ^ 2+m (ti - ^io)-" 1 "' 71 -™ (h - Kt 20 )-^-^~ n \ 

(2.12) 

As indicated, we will use a short-hand notation suppressing the dependence on the vari- 
ables s,ti,t 2 . We now have to choose appropriate contours for the integral representation 
(2.1). For each of the integrals over t\ or t 2 a closed circle, traversed once in the positive 
sense, around the relevant singular point of the integrand is appropriate, since the perti- 
nent exponent is an integer. For the s -integral we need an infinite contour which starts 
somewhere at infinity where the exponential factor of the integral vanishes, surrounds one 
of the singular points in the positive sense, and returns (on a different sheet) to the starting 
point. Assuming that sq is real and positive, then, if 

< arg(z) < iyr (2.13) 
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the starting- and end-point at infinity has the phase tt/2. If we agree that 

arg(s - ksq) = (2.14) 

when s is positive and sufficiently large, the integral of a single term of the infinite series 
(2.4) can be evaluated: 

(nt 2 o+) (rctio+) (kso+) 

z x (2m)~ 3 j> j> J ex.p(z 3 s + z 2 t 1 + zt 2 ) 

/tso+ioo 

X ( S - K S0 ) M(K)+m (tl - Ktl )- Vl - ni (t 2 ~ Ktzoy^-^dsdt^h 2A ' J) 

g— t(k) — 3m+2ni +712 

= eX.p(KSQZ 3 + KtiQZ 2 + Kt 20 z): 



T(-/jl(k) - m)r(i/i + n 1 )Y{v 2 + n 2 ) ' 
where, in the power of z on the right-hand side, already use has been made of (2.5). As 
a consequence, the integral of V(k; qi,q 2 ) yields, if the series is integrated term by term, 
one or the other of the formal solutions (1.2)-(1.3): 

(K,t 2 0+) («tlO + ) ("S0 + ) 

foojiz) :=z x (2ni)- s f f J eM^s + zHi + zh) 

x V(k; q 1 ,q 2 ;s,ti,t 2 )dsdtidt 2 ~ f™*(z) 

in the sector < arg(z) < ^7r, where j = 1 if k = 1 or j = 2 if k = — 1 . Each of 
the solutions defined by the integral representation (2.16 ) has one of the formal solutions 
(1.2) or (1.3) as its asymptotic expansion as z — > oo in the indicated sector. It follows by 
rotation of the contour that the asymptotic expansions are theoretically valid in the larger 
sectors — qTT < arg(z) < |-7r for j = 1 or — ^tt < arg(z) < \n for j = 2, respectively. 

Looking at (2.15) we may see that all the terms for which 3m — 2n\ — n 2 is the same 
yield the same power of z. For the coefficients of the asymptotic expansions (1.2) or (1.3) 
we then have the representation 

a n (K) = ^2 6(K;m,ni,n 2 ), (2-17) 

(m,n 1 ,n 2 )el n 

where 

I n = {(m, ni,n 2 ) : 3m - 2n\ - n 2 = n} C N x N x N , n = 0, 1, (2.18) 
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Because of the properties of the b(m, ni,n 2 ), the sum in (2.18) is finite for each finite n, 
in particular we have 

oo(k) = 6o(k; 0,0,0) = 1. (2.19) 



3. Analytic continuation of the integrand 

Below we have to consider the integral representation (2.1) with t\ and t 2 -contours which 
are simple closed curves surrounding in the positive sense both the finite singular points 
— £io and t\Q or — t 2 o and £20, respectively, and with an s-contour which starts at or near 
—So + ioo, surrounds both the finite singular points — sq and sq once in the positive sense 
and ends at sq + ioo. We therefore need the analytic continuation of the integrand between 
the two singular points along this contour. With appropriate power factors $ included for 
later convenience, the continuation formula reads 
r 1 ,r 2 )V(K;r 1 ,r 2 ) 

00 00 (3 1) 

= ^(- K ; r i' r 2;?i^2)$(-«;?i,?2)V(-K;Qi,Q 2 ) + U(-K;r 1 ,r 2 ), 

<?l=ri q 2 =r 2 

where 

<&(«;ri,r 2 ) = (-2^ r M(K) ^ ri -K (3.2) 

The effect of these power factors is that in (3.1) the total powers with non-integer exponents 
are powers of |[1 — s/(kso)] or of |[1 + s/(kso)], respectively. We may agree that arg(l — 
s/so) = arg(l + s/sq) = when s is on the real axis between — so and sq. Also, as above, 
arg(s) = when s is larger than so- Then, by analytic continuation along the contour 
under consideration, we have 

$(1; n, r 2 ) = (2 So )- M(1) -i ri -^ 2 exp(i7r(^(l) + |n + |r 2 ), (3.3) 

$(-1; n, r 2 ) = (2s )-^- 1) -i ri -^. (3.4) 

Let us rewrite (3.1) using an even more condensed notation, writing q for (<2i,(/ 2 ) in the 
parameter list of the various functions and writing a sum over q in place of a double sum 
over qi and g 2 , etc. The above continuation formula then reads 

r)V( K ; r) = ]T £?(-«; r, g) + r), (3.5) 
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and the second continuation formula 
JJ(k; r) = 3>(-k, r)V(—K, r) 

"EE E ^ r ' 9. p)*(-«; p)^(-«; p) - E ^ r ' ?)^(-«; 9) (3 ' 6) 

follows by the requirement of consistency of (3.5) for k replaced by —k. 

4. Asymptotic expansions for the coefficients in the continuation formula 

We now want to determine the .E-coefficients in the continuation formulas (3.5) , (3.6) by 
means of the asymptotic method of Darboux [16] applied to the variable s . The left-hand 
side of the continuation formula (3.5) is 



oo m m 



(\ - ^r (K)+ ~ 3ri+1 * r2 E E E r ^ + r i+ + r *+ »*) 

m=0ni=0n2=0 

x T(-/j(k) - |n - \r 2 - m)6(«; m, n u n 2 )(-2KS ) m ( 4 ' 1 ) 
x (I " 77^) m (^i " «*io)- , ' 1 - ri - ni (*2 - ^ 20 )-^- r2 - n2 . 

The leading singular term, when (s,£i,£2) — ► (— ft^o, — ^io, —^20)5 on the right-hand side 
is 

00 00 



E E ^(^ K ; r i' r 2;?i^2)r(-^(-K) - |qi - |g 2 ) 

x r(i/i + qi )T(u 2 + q 2 )b{-K- 0, 0, 0) ( 4 - 2 ) 

x (I + -^r { - K)+ i q ^+b^ tl + Kt 10 )-^-^(t 2 + Kt 20 )-^- q2 . 

By means of the binomial theorem in its hypergeometric-series-form 

CO / s 

(i-aO- Q = E^r^ ( 4 - 3 ) 

3=0 J ' 

where 

(a)j = a(a + 1) . . . (a + j — 1) = T(a + j) /T(a) 
means the Pochhammer symbol, we may expand, if \t\ is sufficiently large, 



3=0 J 



(4.4) 



and, if \s/(kso)\ is sufficiently small, 

(| + = f; t±±^z&± (i _ _^ r . (4 . 6) 

m=0 

Then the leading singular term on the right becomes 

oo oo 

Yl Yl r 2 ; qi, 92)r(-//(-K) - fgi - |g 2 ) 

<Jl=i-i 92=^2 

V (^1 + gl)ji (_2^ 10 ) jl (t! ~ Ktw)-"*-*-* 1 

3^0 * ! 

x T(u 2 +q 2 ) Y ( V2+ f )ja (-2 K t 20 y>(t 2 - Kt.20)- 



(4.6) 



32=0 

2„ 1 



x 6( _ K: , ,o) £ - -£_)-. 

^-^ m! z 2ks 

The coefficients of this series should agree asymptotically, as m — > oo, with those of the 
series on the left-hand side, where the power factor in front of the series (4.1), when 
s — > — ksq, tends to 1 and may be omitted. We therefore obtain 



m m 



^ ^ T(z/i + ri + ni)r(z/ 2 + r 2 + n 2 )r(-//(«) - |n - |r 2 - m) 

ni=0 n2=0 

x 6( K ;m,n 1 ,n 2 )(-2^ ) m (ti - /^lo)"" 1- " 1 (* 2 - /^o)"" 2 ""™ 



oo oo 



X] X] E(-K;n, r 2 ;qi, q 2 )T(-fi(-K) - |?i - §g 2 ) 



gi=i"i g2=r2 

oo 



x% + 9l ) ]T (-2^ 10 ) J1 ^i - rtio)-*-*-* (4 ' 7) 



Ji=0 - 71 



x I> 2 +q 2 ) V ( ^ + f )j3 (-2/rf 2 o) Ja (t 2 - Kt 20 y 



■^2-Q2-j2 



32=0 

x 6(-k; 0,0,0) — 



A»(-k) - §<?1 - 592)m 



m! 



which holds asymptotically as m — > oo. On both sides of this asymptotic equation a double 
power series in the same variables appears, so the coefficients of the corresponding terms 
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must be equal. This yields, for each set ri, r 2 , ni, n 2 , 

r(-/i(«) - |n - \r 2 - m)(-2KSo) m b(K; m, m, n 2 ) 

n+nx , 2 +n 2 r(-M-«)-§gi-^ + m) 
« 1. Z. £(-«;ri,r 2 ; 9l ,<fe) ^ 

x (-L(-2k£ 10 )^ f ^- r (-2^t 20 )- 72 ^ &(-«; 0,0,0) 



or, if we introduce new indices of summation and make use of the reflection formula of the 
gamma function, 

- 7r[r(l + m)]~ 1 (2KS ) m b(K; to, n u n 2 ) 

P 7~^o r(l + m) 

v r(l + //(«) + §^i + ^2+m) . 2 ! 

x r(l + m) sm(7r[/x(«) + fl + 3^]) 

x fJ_(-2«t 10 )^ > ) (_L(_ 2/ ^ 20 )^ &(_«; 0,0,0). 



\3l- / ji=ni-pi V-? 2 ' / J2=n 2 -P2 

(4.9) 

Now the left-hand side is independent of ri and r 2 and so is the product of the two ratios 
of gamma functions, when to — > 00, on the right. Therefore £?(— k; ri, r 2 ; ri + p\,r 2 + 
p 2 ) sin(7r[^(«) + |ri + 3 7*2]) must be independent of ri and r 2 too. This proves 
Lemma 3. 

^/ \ sin(7ru(re)) . 

£? -«; n, r 2 ; n + Pl , r 2 + p 2 ) = . 1 ^ JJ i r- ^ -«; 0, 0; Pl ,p 2 ). (4.10) 

sin(7r[/i(«;J + gTi + gr 2 J) 

After Lemma 3 is available, we need to determine the ^-coefficients for T\ = r 2 = 
only, and it is advantageous to introduce the closely related coefficients 

e(-K;p 1 ,p 2 ) = sm(7rfj,(K))E(-K;Q,Q;p 1 ,p 2 ), (4.11) 

which are independent of A. Eq (4.9), with n = r 2 = and the factor b(— «; 0, 0, 0), which 
is equal to 1, omitted, then becomes 

— TV 

— — — — — -(2Ks ) m 6(K; m, rei, n 2 ) 

r(l + f-iyK) + m) 

r(-//(-«) - |pi - |p 2 + to) 

« Z, e(-K;pi,P2) fTTTH ( 4 - 12 ) 

Pl =0p 2 =0 v y 

x f-L(_2/rf 10 ) J ' i> ) ^(-2^20)^) 
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We then may solve (4.12) for the e-coefficient with pi = n\,p 2 = n 2 and obtain the 
following asymptotic formula in terms of a ^-coefficient and the earlier e-coefficients, 



e(— k; ni, n 2 ) ~ — iv(2KSo) m 

r(l + m) 



x 



6(«;m, m, n 2 ) 



r(l + //(«) + m)r(-//(-«) - |ni - |n 2 + m) 

ST ST < \ T (-v(-k) ~ %V\ ~ \P2 + m) 

1^ e (-^P^P2)7T ( f v 2 1 7 

^o P % r(-//(-«) - 5 m - + m) 

(pi,P2)^(ni,n 2 ) 



(4.13) 



x ( -L(-2^ 10 )^ 



Ji=ni-pi 



^7(-2^20) J2 
J2! 



J2=n 2 -P2 



By repeated application of this formula, all the e-coefficients on the right-hand side may 
be eliminated, and this yields the remarkable explicit limit formula 



e(-K;ni,n 2 ) = -ir 



T(l + m) 



x lim 

m^oo r(i + ^(k) + m)r(— «) 



\n 2 + m) 



(4.14) 



m n 2 



x (2^o) m ^ J] K^mJuh) 

Jl=0j2=0 



n 2 ~j2 



(ni-ji)\ (n 2 -j 2 )\ 



The proof of this formula will be given below in Section 7. 

Although approximate asymptotic equations or limit formulas such as (4.13) or (4.14) 
are interesting from a theoretical point of view, we finally need more, namely a detailed 
asymptotic expansion suitable for accurate numerical evaluation. This can be obtained, 
on the basis of Schafke and Schmidt [17], essentially in the same way as above, apart from 
the following two refinements: The power factor in front of (4.1) can no longer be omitted, 
and we have to include a finite number of singular terms rather than the leading one, (4.2), 
alone. The result of this procedure, which will be derived in more detail below, is 

Theorem 1. The E-coefficients, or e-coefficients according to (4.11), in the continuation 
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formula (3.1) or (3.5) are 
e(-K;ni,n 2 ) = —tt 



x 



r(l + m) 



r(l + //(k) + m)r(— fJ,(— k) — gTT-i — 3 77-2 + m 

r(-//(-«) - |gi - ±<7 2 +to 



(2«so) m 6(«; m, ni, 77-2) 



m n 2 



9i=0 g2=0 
(«i,52)^(ni,n 2 ) 

if k k 



r(—fi(—K.) - \n\ - \n 2 + to) 



x 



1 + Z^ Z^ Z^ 7TT ; TT 2 , 1 -H(-K;k,h,l2;qi,q2) + 0(m 



x _(-2«* 10 )* 
Ji' 



Ji+?i+ii=ni V^2 



— (-2^ 20 ) J2 



X 



1 + 2^ iTTT7 TTT2TTTTT ^ /c ' °' ni ' 712 ) + °( m ) 

J (1 + //(-Kj + 3TO1 + 3?l2 - TOjfc 



-1 



(4.15) 



where 



H(-K;k,h,l 2 ;qi,q2) = 



(M«))fc- 



(fe-j)!(l + At(-«) + |?i + 392 



2 , (-2KSoyb(-K;j,hJ 2 )- 



(4.16) 



With a suitable choice of to and iY, Theorem 1 may be used to compute accurate values 
of the e-coefficients, beginning with 



e(— k; 0, 0) = — 7r 



r(l + m) 



r(l + //(k) + m)T(-fi(-K) + to) 



(2^ ) m K«;w,o,o) 



x 



K 



1 + ^ 



(1 + H{-K)) k 

" x (! + //(-«) -m)fc 



#(-«; fc, 0, 0; 0, 0) + 0(m- K_1 ) 



(4.17) 



While the e-coefficients are independent of A, the approximate values computed for any 
finite m do depend on it. Thus A here plays the role of a computational parameter which 
could be adjusted for optimal accuracy. The best choice from this point of view , in 
particular when L is not small, is A = L or A = —L, according to the discussion in Section 
7 below. 



5. Multiplicative solutions and the Floquet exponent 
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We now want to construct a linear combination of the solutions at infinity, 

f p (z) = a/ooi(z) + (5foo2{z), (5.1) 

which is a multiplicative solution such that, after analytic continuation along a sufficiently 
large circle around the origin traversed once in the negative sense, this solution remains 
the same apart from multiplication by a constant factor, that is 

f p (e- 2 ^z)=pf p (z). (5.2) 

This solution is proportional to one of the Floquet solutions introduced above in the 
introduction. 

We use the integral representation with contours which surround both the finite singular 
points as introduced above in Section 3. This integral is equal to the sum of the two 
integrals with contours, also considered above in Section 2, which surround only one of 
these singular points. If the contour is kept fixed, then the circle in the z-plane traversed 
once in the negative sense corresponds to a circle, around the origin and with radius greater 
than so, in the s-plane traversed three times in the positive sense, since then z s s in the 
exponential part of the integrand does not change. 

Let us consider the integral representation with V^(1;0) in the integrand and with the 
above phase conventions. From (3.5) we have 

$(1; 0)^(1; 0) = 0, q)V(-l; q) + U(-l; 0). (5.3) 

q=0 

The integral then yields a function 

f\z) := $(1; O)/^) + 0, q)foo 2 (z), (5.4) 

q=0 

where the two terms come from the two singular points, which contribute via the left- 
or right-hand side of (5.3), respectively. Let us now consider, in the s-plane, a simple 
closed loop, homotopic to the circle mentioned above, consisting of small circles around 
the singular points sq and — sq and straight line segments along the real axis between these 
singular points. Let us start with the continuation formula (5.3) in a neighbourhood of 
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the origin, where the left- as well as the right-hand side of (5.3) are valid, and see what 
happens when we follow the loop in the positive sense. We shall always refer to the above 
phase conventions and display any additional phases explicitly. Traversing the circular 
part around sq multiplies the function 1/(1; 0) by the phase factor exp(27ri / u(l)), so that 
we then have 

exp(27ri^(l))$(l;0)y(l;0) 

= exp(27ri^(l)){^ E(-l; 0, q)$(-l; q)V(-l; q) + U(-l; 0)}, (5 ' 5) 

q = 

where the right-hand side is the analytic continuation of the left by means of (5.3). Next 
we have to traverse the circular part around — so, which multiplies the function V(—l;q) 
by the phase factor exp(27ri[ / u(— 1) + 99(g)]), where 

¥>(?) = h 1 + l^ 2 ' ( 5 - 6 ) 

so that we have on the straight line segment from — sq to sq , after the loop has been 
traversed once, 

exp(27ri[Ml) + M-1)D E E ^ °> 9) exp(27ri^(g))$(-l; q)V(-l; q) 

q=0 

+ exp(27ri/x(l))t/(-l;0) 
= exp(27ri^(l)){$(l; 0)^(1; 0) + ^ {exp(27ri[^(-l) + <p(q)]) - 1}U(1; q) 

(5.7) 

+ EE E (-^ °' 9){exp(27ri[^(-l) + <p(q)}) - 1}E(1; q,p)Q(l;p)V(l; P )}, 

q=0 p=q 

where again the right-hand side is the analytic continuation of the left. It is evident that 
following the loop further leads to increasingly lengthier and more complicated formulas, 
not suitable for being fully displayed in this paper. We therefore want to stop here for a 
moment and consider the integral representation with the integrand obtained after the loop 
in the s-plane has been traversed only once (rather than three times, as finally needed). A 
representative example of the terms in (5.7) then is 

E5]£;(-l;0,g){e X p(27ri[^(-l) + ^)])-l}£;(l;g,p)$(l;p)y(l;p). (5.8) 
q=o p=q 
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The integral of V(l;p) yields fooi(z), independent of p, so that the sum over p can now 
be performed and , because of Lemma 3, the multiple sum reduces to a product of single 
sums. After integration we therefore obtain for (5.8) 

{2isin(7r/i(-l)) exp(i7r/i(-l)) ^£(-1; 0, g)(2s )~^ (9) exp(2niip(q)) 

(5.9) 

x ^E(l;0,p)(2 S o)-^ ) exp(i7r V? (p))}$(l;0)/ ool (z). 

p=0 

It is now convenient to introduce the phase factor 

?? :=exp(±7ri) = §(l + h/3) (5.10) 

which satisfies 

77 6 = 1, 1 + V 2 + V 4 = 0. (5.11) 

Also, for the sums occurring here and below we may introduce the "Stokes multipliers" 
<r n (l) :=sin(7r^(l))^£;(-l;0, 9 )(2 S o)-^ (,?) exp(2n7ri^( 9 )), 



(5.12) 



(5.13) 



q=0 

a n (-l) := sin(7r/i(-l)) J]E(1; 0, q)(2s )-^ exp((2n + l)Tri^(g)), 

<?=o 

which become 

a (l) = 5 (1) + (2 So )- 1/3 Si(l) + (2 So )- 2/3 ^ 2 (l), 

ao(-l) = So(-l) + v(^o)- 1/3 S 1 (-l) + 77 2 (2 So )- 2/3 ^2(-l), 

ffi(l) = 5 (1) +r7 2 (2 So )- 1/3 ^(l) + r7 4 (2 So )- 2/3 ^(l), 

(n(-l) = 5 (-l) + r7 3 (2 S o)- 1/3 ^i(-l) + (2s )- 2/3 S 2 (-l), 

a 2 (l) = 5 (1) +r 7 4 (2 So )- 1/3 S'i(l) + 77 2 (2 So ) " 2/3 ^2 (1) , 

<r 2 (-l) = 5 (-l) + r 7 5 (2so)- 1/3 S' 1 (-l) + 7 ? 4 (2 So )- 2/3 ^(-l) 
in terms of the real partial sums 

oo 

S (k) :=^2(2ks )~ 1 ^2 e(-/c;ni,n 2 ) = e(-«; 0,0) + . . . , 

^=0 (ni ,rt 2 )e J3i 

oo 

Si(k) :=^2(2ks )~ 1 e(-K;ni,n 2 ) = e(-«;0,l) + ..., (5.14) 

z =0 (ni,n 2 )£ J3i+i 

00 

5 2 (k) :=^2(2ks )~ 1 Yl e(-K;m,n 2 ) = e(-K;l,0) + e(-K;0,2) + ..., 

'=0 (ni,n 2 )e J3Z+2 
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where 

Ji = {(n 1 ,n 2 ) : 2m + n 2 = 1} C N x N , / = 0, 1,.... (5.15) 
The integral of the representative term (5.9) then becomes 

2iexp(i7r / u(— 1)) 



sin(7r/i(l)) 

In total, we have 



<7i(l)"o(-l)*(l;0)/ooi(*). (5.16) 



sm^7T//^l j j 

+ exp(27ri^(l){l + 2ie s ^^~ 1)) a 1 (l)a (-l)}$(l; 0)/ooi(s), 



(5.17) 



where 



/ T (z) = $(1; 0)/ooi(z) + — 7^(T)) ^(1)$(-1; 0)/oo2 (z). (5.18) 



So far we have traversed the s-loop once and obtained on it the analytic continuation of 
the integrand, but we have to traverse it three times. This yields 



f l {er im z) = Su$(l; Oj/ooiCz) + 5 12 $(-1; 0)/ oo2 (^) (5.19) 

with lengthy expressions for Su and S±2- 

In a similar way, the integral representation with V{— 1; 0) yields 

/ n (^) := ]T £(1; 0, g)/ooi(*) + *(-l; 0)/oo 2 (*) (5.20) 

q=0 

or 

(*) = ■ 7 \ ^ ^(-1)^(1; O)fooi(z) + $(-1; 0)/oo 2 (z) (5.21) 
sin(7r/i( — lJJ 

and leads, after the loop has been traversed once, to 

exp(27riA/3)/ n (e- 2 -/^) = {2i exp(7ri^(-l))a (-l) + g^^k i(-l) 

sin (717^— 1 j j 

_ 4 e*p W Ml) +,(-!)]) 
sm(7T/x(— 1)) 

+ ex P (27ri^(-l)){l+^^^(7 1 (-l) ( 7i(l)}*(-l;0)/ O o2(z) 
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and, after three times, 

f u (e- 2 * [ z) = S 21 $(l; O)/^) + S 22 <&(-1; 0)/ oo2 (z), (5.23) 

where again the expressions for S21 and S' 22 are too lengthy to be displayed here. 
What we really want to obtain are the circuit relations for /^i and /002, that is 

/ooi(e- 27 ^) = Tu/oci^) +212/002(2), 

/oo2(e" 27r1 ^) = T 2 i/ooi(z) +T 22 / oo2 (2), 

where 

r n = - _4»a_ S21 ,/[i -o(-D-o(i) , 



(5.24) 



sin(7r / u(l)) sin(7r / u(— 1)) sin(7r / u(l)) 

T 12 = [$(-1; 0)/*(l; 0)] [5 12 - S 22 ]/[l -o(-l)-o(l) 



sin(7r / u(l)) sin(/T / u(— 1)) sin(7r/x(l)) 

T 21 = [*(!; 0)/*(-l; 0)] [5 21 - . ^W n]/[1 ^(-l)ao(l) 



(5.25) 



sin(7r / u(— 1)) sin(7r / u(— 1)) sin(7r/i(l)) 



r 22 = [& - , a ; ( -V s 12 ]/[i - -"{-Doom 



sin(7r/i(— 1)) sin(7r / u(— 1)) sin(7r / u(l)) 

It turns out that each numerator contains a common factor which compensates the de- 
nominator. The result then is 

Theorem 2. The coefficients in the circuit relations (5.24) are 
T n = exp(27rir(l)) 

+ 4exp(|7rir(-l))(To(-l)(To(l)+4exp(|7rir(l))(Ti(l)(To(-l) +4a 2 (l)(To(-l) 
+ 4 exp (|7rir(l))(7 2 (l)(7 1 (-l) + 4 exp (§7nr(l))c7 2 (-l)c7 (l) + 4a 1 (-l)a (l) 
+ 16exp(|7rir(l))a 2 (-l)(T 2 (l)(7i(-l)(To(l) + 16 exp(§7nr(l))c7 2 (-l)c7i(l)c7o(-l)ao(l) 
+ 16exp(§7rir(l))(7 2 (l)(7 1 (-l)a 1 (l)(7o(-l) + 16 exp(§7riT(-l))<7 2 (-l)a 2 (l)<7 (-l)<7o(l) 
+ 16exp(|7rir(-l))(7 1 (-l)(7 1 (l)(7o(-l)(7o(l) 
+ 64a 2 (-l)a 2 (l)(7i(-l)ai(l)(7o(-l)(7o(l), 

T 12 = (2 So ) (1/3)[T(1) - r( - 1)] {2iexp(|7rir(-l))(7 2 (l) 
+ 2iexp(|7rir(l))(Ji(l) + 2iexp(27rir(-l))a (l) 

+ 8iexp(|7rir(-l))(7 2 (-l)a 2 (l)(7o(l) + 8i(7 2 (-l)(7i(l)(7o(l) (5.266) 
+ 8iexp(|7rir(-l))(7 1 (-l)a 1 (l)(To(l) + 8ia 2 (l)a 1 (-l)(7 1 (l) 
+ 32iexp(|7rir(-l))(7 2 (-l)(7 2 (l)(7 1 (-l)(7 1 (l)(7o(l)}, 
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(5.26a) 



T 21 = (2 So ) (1/3)[r( - 1) - T(1)] {-2iexp(|7rir(-l))a (-l) 

- 2iexp(|7rir(l))a 2 (-l) - 2kr 1 (-l) 
-8iexp(|7rir(l))(7 2 (-l)(7 2 (l)(7 1 (-l) 

- 8iexp(§7rir(-l))(T 2 (-l)(7 2 (l)(To(-l) (5.26c) 
-8iexp(§7rir(-l))(7 1 (-l)(7 1 (l)(7o(-l) 
-8iexp(§7rir(l))(7 2 (-l)(7 1 (l)(7o(-l) 
-32ia 2 (-l)(7 2 (l)a 1 (-l)(7 1 (l)ao(-l)}, 

T 22 = exp(27rir(-l)) 

+ 4exp(|7rir(-l))(T 2 (-l)(T 2 (l) +4exp(|7rir(-l))ai(-l)ai(l) +4a 2 (-l)ai(l) (5.26d) 
+ 16exp(§7rir(-l))(7 2 (-l)(7 2 (l)(7 1 (-l)(7 1 (l). 
It turns out that 

TnT 22 -T 12 T 21 = l. (5.27) 

We are looking for a multiplicative solution such that (5.2) holds. It then follows from 
(5.1), 5.2), (5.24) that 

(Tn - p)ct + T 21 p = 0, 



T 12 a + {T 22 -p)P = 



(5.28) 



and, as a consequence, that 



(T 11 -p)(T 22 -p) = T 12 T 21 (5.29) 



or 



p 2 -(T n +T 22 )p + l = 0, (5.30) 
where (5.27) has been used. The roots pi, p 2 of this equation satisfy 

Pi+P2 = T n +T 22 , (5.31) 

and they may be represented in terms of one (not necessarily real) parameter u as 



p\ = exp(— 27rio;), p 2 = exp(27ria;). 
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(5.32) 



Then we have 

Pi +P2 = 2cos(27ro;) = Tn + T 22 , (5.33) 

and the final result is 

Theorem 3. The characteristic exponent u of the multiplicative solutions is given by 

cos(2ttw) = cos(2ttt(1)) + X, (5.34) 

where 

X = 2exp (f 7riT(-l))<7 (-l)<7 (l) + 2exp (§7rir(l))<7i(l)<7o(-l) + 2<t 2 (1)<7 (-1) 
+ 2exp(|7rir(l))(7 2 (l)(7i(-l) + 2exp(|7rir(l))a 2 (-l)(7o(l) + 2(7i(-l)(7o(l) 
+ 2exp(|7rir(-l))(7 2 (-l)(7 2 (l) + 2exp(|7rir(-l))ai(-l)(7i(l) + 2a 2 (-l)(7i(l) 
+ 8exp(§7rir(-l))(7 2 (-l)(7 2 (l)(7i(-l)(7i(l) 

+ 8exp(|7rir(l))(7 2 (-l)(7 2 (l)ai(-l)(7o(l) + 8exp(|7rir(l))(T 2 (-l)(7i(l)(7o(-l)(To(l) 
+ 8exp(|7rir(l))(7 2 (l)(7i(-l)(7i(l)(7o(-l) + 8 exp(|7rir(-l))<7 2 (-l)(T 2 (l)<7o(-l)(To(l) 
+ 8exp(|7rir(-l))ai(-l)(7i(l)(7o(-l)(7o(l) 

+ 32(7 2 (-l)(7 2 (l)(7l(-l)(7l(l)(7o(-l)(7o(l). 

(5.35) 

For each of the roots p\ and p 2 we may determine the ratio of a and (3 from the upper 
or lower equation of (5.28). Choosing in each case a convenient normalization, we get the 
desired multiplicative solutions 

f pl (z) = {T 22 -p 1 )f ool {z)-T 12 f oo2 {z), 

(5.36) 

f P 2(z) = (T 22 -p2)/ool(*) -T 12 foo 2 (z). 

Here the lower equation of (5.28) has been used in both cases so that T 11: which consists 
of a considerably longer expression than T 22 , does not appear. 

Whenever D m = for all the m = 1,2, . . . , 6, the origin is a regular singular point 
of the differential equation with exponents uj = L, —L. We are not able, however, to see 
analytically that then the lengthy expression (5.35) reduces to X = cos(2-7rL) — cos(27rr(l)), 
but this is confirmed in examples of numerical computations, as expected. 
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6. Asymptotic behaviour of the late coefficients of the formal power series 
solutions 

According to (4.12), the leading terms of the asymptotic behaviour of the ^-coefficients 
for large to are given by 

6(k; to, m,n 2 )« 2kso Wi v — - 

7T r(i + m) 

x ^ ^ e(-K;pi,p2)r(-A»(-K) - fpi - §P2 + m) (6.1) 

Pi=0 p 2 =0 

(ni-pi)! (W2-P2)! 
This result may be used to discuss the asymptotic behaviour of the late a-coefncients of 

the formal solutions. Writing the decomposition (2.17) of the a-coefncients in terms of the 

6-coefncients separately for each of three consecutive indices n = 3iV, 3N + 1, 3N + 2, we 

have 

2N 

a 3 N(K)=J2 b(K;N + l, ni ,n 2 ) = 6(k; AT, 0, 0) + . . . , (6.2) 

1=0 (n 1 ,n 2 )eJ3l 

2N+1 

a 3N+1 (K)=J2 Yl b(K;N + l,n 1 ,n 2 ) = b(K;N + l,l,0) + b(K;N+l,0,2) + ..., 

1 = 1 (n 1 ,n 2 )eJ 3 l-l 

(6.3) 

2N+2 

a 3 N+2(K)=J2 Yl b(K;N + l, ni ,n 2 ) = b( K ; N + 1,0,1) + ..., (6.4) 

1 = 1 (n 1 ,n 2 )eJ 3 i-2 

where 

Ji ={(ni, n 2 ) : 2n x + n 2 = 1} C N x N , / = 0, 1,.... (6.5) 

Inserting (6.1) and omitting terms of relative order iV -1 , we get from (6.2) 

r \ n ( \ V"" ^ r M-^-n)-i Pl -lp 2 + l + N) 
a 3 N{K) ~ C 3N (k) ^ ^ 2^ 2^ e(-/c;pi,p2) rf-//f-«) + AO 



x (2^o)" Z 7 ^(^xo)" 1 -^, ^— (-2^ 2 o)" 2 - p2 , 

(ni-pi)! {n 2 -p 2 y. 



where 

r /x 1 r(l + h(k) + N)T(-ijl(-k) + N) N 
« --r(/x(«) - /i(-/c) + N)(2ks )~ N . 

7T 
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(6.6) 



(6.7) 



Let us now discuss (6.6) in detail: Here only such values of n\ and n 2 occur for which 

§ni + \n 2 = I. 

If again terms of relative order iV -1 are omitted, the ratio of the gamma functions is equal 
to 

N -(2/3) Pl -(l/3) P2 +l _ jy(2/3)(ni-pi)+(l/3)(n 2 -P2)_ 

The asymptotic behaviour of (6.6) then becomes 

27V ni n 2 

1=0 (n 1 ,n 2 )eJ 3 l Pi=0p2=0 

x [iV 2 / 3 (2^ )- 2/3 (-2^io)] ni - pi 7 ^- ? [iV 1 /3 ( 2 KSo )-V3 ( _ 2 ^ 2o)] n 2 - P2 _ 

(ni-pi)! (n 2 -p 2 )! 

(6.8) 

Here the integer power of 2ksq has been splitted in two factors with fractional powers, the 
meaning of which is given in terms of the phase factor r\ of (5.10) by 

(2K8) -i_ f(2 So )- (2/3) -(2 S o)- (1/3) - if « = 1 (6Q) 

With n/3 in place of AT, the terms in (6.8) look like the terms of the expansion of the 
exponential function 

EX (n) : = exp(( ±n) 2 / 3 (2/^ )~ 2/3 (-2^io) + (§7i) 1 / 3 (2ks )~ 1/3 (-2k*2o)) 

(6.10) 

= 1 + (in) 2 / 3 (2^ )- 2 / 3 (-2^ 10 ) + (in) 1 / 3 (2^ )- 1 /3 ( _2^ 20 ) + . . . , 
but because of the three possible values of a third root there are two other such functions, 
EX 1 (n) : = exp(r ? 4 (in) 2 / 3 (2^ )- 2 / 3 (-2^ 10 ) + V 2 ^n) 1 / 3 (2 K s )- 1 / 3 (-2 K t 20 )) 

= 1 + ?7 4 (in) 2 / 3 (2^o)- 2/3 (-2^io) + v 2 ^n) 1 / 3 (2KS )- 1/3 (-2Kt 20 ) + . . . , 

(6.11) 

EX 2 (n) : = exp(? ? 2 (in) 2 / 3 (2^o)- 2/3 (-2^io) + v 4 ^n) 1 / 3 (2KS )- 1/3 (-2Kt 20 )) 

= 1 + r ? 2 (in) 2 / 3 (2^ )- 2/3 (-2^io) + r ] \\n) 1 / 3 {2KS )- ^ 3 {-2nt 2Q ) + . . . . 

(6.12) 

Therefore (6.8) asymptotically shows exponential behaviour given by a certain linear com- 
bination of these three exponential functions. They can be identified by the constant term 
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and the two linear terms of their expansion shown above. It is convenient, and easy because 
of the properties (5.11) of rj, to introduce three functions which are linear combinations 
of the exponential functions such that only one of the identifying terms is different from 
zero, 

L (n) := \[EX (n) + EX 1 {n) + EX 2 (n)] = 1 + . . . , 

Li(n) := \{EX (n) + rfEX^n) + V 4 EX 2 (n)} = (±n) 2 / 3 (2/^ )~ 2/3 (-2^io) + • • • , 

L 2 {n) := \[EX {n) + r ] 4 EX 1 (n) + rfEX^n)] = (^n) 1 / 3 (2KS )- 1 / 3 (-2Kt 20 ) + .... 

(6.13) 

We can now determine the asymptotic behaviour of (6.8) by looking at the contributions 
from {pup 2 ) = {ni,n 2 ) or {pup 2 ) = (ni - l,n 2 ) or (pup 2 ) = (n u n 2 - 1), respectively. 
This yields 

a 3 iv(«) ~ C 3N (k) {s (k)L (3N) + (2ks )- 1/3 S 1 (k)L 1 (3N) + (2^ )" 2/3 S 2 (k)L 2 (3N)} . 

(6.14) 

In a similar way we may obtain from (6.3) and (6.4), respectively, 
fl3JV+l(«) ~ C 3N+1 (k) 

x {(2ks )- 2/3 S 2 (k)L (3N + 1) + S (ft)Li(3iV + 1) + (2^ )- 1/3 ^i(k)L 2 (3A^ + 1)} , 

(6.15) 

x {(2^o)- 1/3 5i(k)L (3AT + 2) + (2ks )~ 2/3 S 2 (k)L 1 (3N + 2) + S (k)L 2 (3N + 2)} . 

(6.16) 

If we now switch back to a representation in terms of the exponential functions (6.10)- 
(6.12), Stokes multipliers, according to (5.13) above, appear as their factors, 

C1 3N (k) S3 ^C 3N (k) 

x {cto(k)EXo(3N) + cti(k)EXi(3N) + ct 2 (k)EX 2 (3N)} , 
ci 3N+1 (k) « \C 3N+1 (k) 

(6.17) 

x {a (K)EX (3N + 1) + 77 2 <7i(/c)EXi(3W + 1) + 7i 4 <t 2 (k)EX 2 (SN + 1)} , 

a3Af+2(«) ~ §CW+2(«) 

x {(7 (k)£;Xo(3^ + 2) + 7/Vi(K)£Xi(3iV + 2) + t/V 2 («)£7X 2 (3JV + 2)} . 
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These three asymptotic equations can be combined to give 

Theorem 4. The leading terms of the asymptotic exponential behaviour, as n — > oo, of 
the coefficients of the formal solutions are given by 

a n («) « -i7r- 1 (2^ )" n/3 r(i[r(ft) - r(-«) + n]) 

(6.18) 

x [a ( K )EX (n) +rj 2n a 1 ( K )EX 1 (n) + n 4n a 2 (K)EX 2 (n)]. 

This is a concrete example, with all the quantities determined explicitly, of the structural 
results obtained by Immink [6] , and it is interesting also in the context of related work by 
other authors [5, 15]. 

7. Postponed proofs 

7. 1 Proof of Lemma 1 
Let 

f(z) = z x u(z) (7.1) 

and multiply the differential equation for u(z) by z 6 in order to remove all the negative 
powers of z. We then are concerned with the differential equation 



8„." i /o\ i i\„7„./ 

(7.2) 



L z u(z) :=z*u" + (2A + l)z'u 

6 6 

- [J2 D m z 6 ~ m + (L 2 - X 2 )z 6 + B m z 6+m ]u(z) = 0. 



m=l m=l 



We are looking for a solution of this differential equation in the form of an integral repre- 
sentation 

u(z) = III K{z;s,t\,t2)v{s,t\,t2)&s&t\&t2 (7.3) 

with the kernel 

K = K(z; s, t u t 2 ) = exp(z 3 s + z 2 t x + zt 2 ). (7.4) 

If we perform the differentiations with respect to z under the integrals, the differential 
equation becomes 

L z u = JJJ vL z Kdsdt 1 dt 2 = 0, (7.5) 
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where 

\-L z K = (9s 2 - B G )z 12 + (12sh - B 5 )z 11 + (6st 2 + At\ - B 4 )z 10 



+ ([6A + 9]s + At x t 2 - B 3 )z 9 + ([4A + 4]h + t 2 2 - B 2 )z 

6 

+ ([2A + l]t 2 - B x )z 7 + (A 2 - L 2 )z 6 - °mz 



(7.6) 

6 

6 — m 



m=l 

It is advisable to rewrite this in terms of the quantities p\, p 2 , ps, according to (1.5), 
which determine the exponential factors of the formal solutions (1.2)-(1.3), or in terms of 
the related quantities so, tio, ^20 according to (2.3). Then (7.6) becomes 

1 



—L Z K = 9(s 2 - s 2 )z 12 + 12(s*i - sot^z 11 + (6[st 2 - s t 20 ] + A[t\ - t 2 l0 ])z 10 
A 

+ ([6A + 9]s + 4[*i* 2 - t w t 20 ] + [4t 10 t 20 - B 3 ])z 9 
+ ([4A + 4]*! + [t 2 - t 2 ] + [t 2 - B 2 ])z 8 

6 

+ ([2A + l]t 2 - B^z 7 + (A 2 - L 2 )z 6 - D m z 6 ~ m . 



(7.7) 



m=l 

By repeated partial integrations of the exponential function with respect to s or t\ or t 2 , 
respectively, it is possible to get rid of all the powers of z. This lengthy task can more 
conveniently be performed in a formal way as follows: Let us find a partial differential 
expression M = M s tl t2 with respect to s, ti, t 2 , independent of z, such that 

L z K = M sMM K. (7.8) 

The powers of z correspond to partial derivatives with respect to s or t\ or t 2 , and there 
seems to be some ambiguity as to the choice of the partial derivatives which yield the 
same total power of z. This ambiguity is resolved by the following reasonning: In order 
to get the order of each derivative as small as possible, we would prefer derivatives with 
respect to s, which account for z 3 , with highest priority, next ti, which accounts for z 2 , 
last t 2 , which accounts only for z. Conflicting with this policy, however, are some other 
requirements which ensure that we get a differential expression appropriate for our purpose. 
So the singularities should be at the right places, which are related to the coefficients of the 
exponential factor of the formal solutions we want to represent. In particular, bilinear or 
quadratic factors, such as t\ for instance, should occur in the form t\ —t 2 . For this reason 
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we have in (7.7) already added and subtracted the term At\ . It is then necessary to treat 
these two terms differently, so that any power is always multiplied by the corresponding 
derivative. In the example just mentioned this means that z 10 has to be translated into 
d 4 /ds 3 dt 2 for one term, but into d /ds 2 dt 2 for the other. The unique result now is 

M sMM = 9(s 2 - s 2 )d 4 /ds 4 + 12( S *! - s t 10 )d 4 /ds 3 dt 1 

+ 6(st 2 - s t 20 )d 4 /ds 3 dt 2 + 4(tj - tj )d 4 /ds 2 dtj 

+ ([6A + 9]s + [4t 10 t 20 - B 3 ])d 3 /ds 3 + 4[ht 2 - t 10 t 20 ]d 4 /ds 2 dt 1 dt 2 

(7.9) 

+ ([4A + 4]*! + [t 2 20 - B 2 ])d 3 /ds 2 dt 1 + [t 2 - t 2 20 ]d 4 /ds 2 dt 2 

+ ([2A + l]t 2 - B 1 )d 3 /ds 2 dt 2 + (A 2 - L 2 )d 2 /ds 2 
- D 1 d 2 /dsdt 1 - D 2 d 2 /dsdt 2 - D 3 d/ds - D 4 d/dh - D 5 d/dt 2 - As- 
Next let us introduce the adjoint differential expression M = M S)tl)t2 defined, with any 
sufficiently differentiable function v = v(s,ti,t 2 ), by 

M sMM v = (d 4 /ds 4 ){9(s 2 - s 2 )v} + (d 4 /ds 3 dt 1 ){12(st 1 - s t 10 )v} 
+ (d 4 /ds 3 dt 2 ){6(st 2 - s t 20 )v} + (d 4 /ds 2 dt 2 ){4(t 2 - t 2 )v} 

- (d 3 /ds 3 ){([6\ + 9}s + [At 10 t 20 - B 3 ])v} + (d 4 /ds 2 dt 1 dt 2 ){A(t 1 t 2 - t 10 t 20 )v} 

- (d 3 /ds 2 dt 1 ){([4\ + 4] tl + [t 2 - B 2 ])v} + (d 4 /ds 2 dt 2 ){(t 2 - t 2 )v} 

- (d 3 /ds 2 dt 2 ){([2\ + l]t 2 - B x )v} + (A 2 - L 2 ){d 2 v/ds 2 ) 

- D^v/dsdh) - D 2 (d 2 v/dsdt 2 ) + D 3 (dv/ds) + D^dv/dh) + D 5 (dv/dt 2 ) - D 6 v. 

(7.10) 

The difference as compared with M is that the factors in front of each derivative have 
here to be differentiated too and that all the terms of odd order change their sign. The 
usefulness of the adjoint expression lies in the formula, known in the one-variable case as 
the identity of Lagrange, 

vMK - KMv = RHS, (7.11) 

where the right-hand side RHS is a lengthy expression, bilinear in v and K or their partial 
derivatives, which may be so arranged that it consists of a sum of terms of which each is 
a total derivative with respect to one of the variables. For the term with the factor —D\, 
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for instance, this reads 

d 2 K d 2 v ! d OK _ dv l d OK _ dv 

v lhbT 1 ~ d^dh~ 1 d~s [v dh~ at? + 2 dT 1 [v ~dF ~ ~b~s 

By means of (7.5) and (7.8), the equation to be satisfied now is 



L z u = / / / vMs^Kdsd^dh 



-III 
-III 



KM sMM vdsdt 1 &t 2 + J J J RHSdsdt^h = 0. 



(7.13) 



The first term in the second line can be made to vanish if v is required to be a solution of 
the partial differential equation 

M a , tllt2 i;(*»*i»*2) = 0, (7.14) 

and the second term by a suitable choice of the contours of integration, for it is a sum of 
semi-integrated terms, each involving the difference of the values of the integrand at the 
termini of the contour of one variable and only two remaining integrals with respect to 
the other two variables. The partial differential equation (2.2) is the same as (7.14), after 
the derivatives of the products in (7.10) have been resolved. This completes the proof of 
Lemma 1. 

7. 2 Proof of Lemma 2 
In terms of the shifted variables 

S = s - ks , T 1 =t 1 - ntio, T 2 = t 2 - nt 2 o, (7-15) 

the differential equation (2.2) reads 

9S(S + 2Ks )(d 4 v/dS 4 ) + 12(STi + Kt 10 S + Ks T l )(d 4 v/dS 3 dT 1 ) 

+ 6(ST 2 + Kt 20 S + KS T 2 )(d 4 v/dS 3 dT 2 ) + 4Ti(Ti + 2 K t 10 )(d 4 v/dS 2 dT?) 

+ T 2 (T 2 + 2 K t 20 )(d 4 v/dS 2 dT 2 2 ) + A{T X T 2 + nt^Tr + K t 10 T 2 )(d 4 v/dS 2 dT 1 dT 2 ) 

+ ([81 - 6\}ks + B 3 - 4t 10 t 2 o)(d s v/dS 3 ) + (81 - 6X)S(d s v/dS 3 ) 

(7.16) 

+ ([52 - 4A]k£i +B 2 - t 2 )(d s v/dS 2 dT 1 ) + (52 - 4X)T 1 (d s v/dS 2 dT 1 ) 
+ ([25 - 2\}Kt 20 + B 1 )(d 3 v/dS 2 dT 2 ) + (25 - 2X)T 2 (d 3 v/dS 2 dT 2 ) 
+ ([A - 12] 2 - L 2 ){d 2 v/dS 2 ) - D 1 (d 2 v/dS0T 1 ) - D 2 (d 2 v/dS0T 2 ) 
+ D 3 (dv/dS) + D^dv/dT^ + D 5 (dv/dT 2 ) - D 6 v = 0. 
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Inserting a power series solution 

v = E E E a(m,n 1 ,n 2 )S^ +m Tr i - ni Tr 2 - n2 , (7.17) 

m=0 ni=0 n2=0 

we obtain for the coefficients the recurrence relation 
6ks (/j + m)(fi + m — 1)(/j + m — 2)[3(/z + m) — 2(z/i + m) - (z/ 2 + n 2 ) - A - t(k) + 6] 
x a(m, ni, 77,2) 

+ (/j + m- l)(/j + m- 2){[3(/j + m) - 2(i/i +m) - (v 2 + n 2 ) - A + 3] 2 - L 2 } 
x a(m — 1, ni, 77-2) 

+ + m- 1)0 + m- 2)(i/i +ni - l){4K*io[-3(At + m) + + m) + (u 2 +n 2 ) 
+ A - 4] + t| - S 2 }a(m - 1, m - 1, n 2 ) 

+ (/j + m - 1)(/j + m - 2)(z/ 2 + ?7 2 - l){2Kt 20 [-3(fi + m) + 2{v\ + m) + (u 2 + n 2 ) 

+ A - |] + £?i}a(m - 1, m, n 2 - 1) 

+ -Di(/-t + m — 2)(ui + ?7 1 — l)a(m — 2, 771 — 1, n 2 ) 

+ D 2 (n + 777 - 2)(z/ 2 + ?7 2 - l)a(m - 2,771,772 - 1) 

+ Ds(n + m — 2)a(m — 2, 771, 772) — D^{y\ + n\ — l)a(m — 3, ni — 1, 772) 
- D^{i> 2 + n 2 - l)a(m - 3, 771, n 2 - 1) - D 6 a(m - 3, 771, n 2 ) = 0, 

(7-18) 

valid for m > 0, 771 > 0, n 2 > provided that we agree that all the a-coefficients are equal 
to zero if any of the indices 777,771,772 is less than zero. Assuming that a(0, 0, 0) 7^ 0, we 
get from the equation for m = n\ = n 2 = the indicial equation 

-1)0- 2) (3/j - 2z/i - v 2 - A - r («) + 6) = 0, (7.19) 

with t(k) according to (1.6). Possible values of the exponent \i are therefore ji = 0, 1, 2, or 

/j = § [2^i + v 2 + A + r(/c)] - 2. (7.20) 

In order to avoid complications, we may assume that the last possibility does not yield an 
integer value. This can always be guarantied be a suitable choice of the still disposable 
parameter A. Inspection of the recurrence relation then shows that each of the possible 
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values of [i leads to a solution of the partial differential equation. We here need not 
further consider the three solutions which are regular with respect to S at S = 0, but we 
continue to discuss the singular one with the exponent (7.20), writing for the associated 
coefficients m, ni, 712) rather than a(m, ni, 712). Simplifying by means of (7.20) and 
introducing the 6-coefficients according to (2.6) , we get (2.7). The coefficients for which 
3m — 2ni — n 2 = are constants of integration and may be chosen to be zero. This 
completes the proof of Lemma 2. 

7. 3 Proof of the limit formula 

We have to verify that the limit formula (4.14) satisfies (4.13 ). Substituting it for the 
e-coefficients in (4.13) and interchanging the summations, we have to evaluate sums such 
as 



^{p-j)\{n-p)\ (n-j)\^ Q { } \ q J \ if < j < n - 1, n > 



(7.21) 

Therefore only one term survives on the right-hand side, which then becomes equal to the 
left. 

7.4 Proof of Theorem 1 

In order to prove Theorem 1, we first multiply the continuation formula (3.5), with 
n = r 2 = 0, by 

u j=o J u 

where the left-hand side is used on the left and the right-hand side on the right of the 

continuation formula. Then the left-hand side is the same as above (4.1) with the power 

factor in front of the series removed, but the right becomes, after multiplication of the two 

power series, 

°; ?i» g 2 )r(-M-«) - |?i - b*) 

oo k k 

X J2J2J2 r ( z/1 + 9i + /l ) r ( z/2 + ^ + k > h, h; qu 92) (7 - 23) 

k=o h=o i 2 =o 



x a + _^_)M-«)+f 9i+i*.+fc( tl + K t 1Q )-^-^(t 2 + Kt 20 )- iy2 - q2 - l \ 

ZKSo 
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with H defined in (4.16). Proceeding as above, and making use of the formula 

(a - k) m = (a) m n {1 ~ a) \ , (7.24) 
(1 - a - m) k 



we obtain 

— TV 



■(2KSo) m b(K; m, ni, riz) 



sin(7r/x(«;)) T(l + + m) 

ST ST J? f n n , T(-^(-/t) - jgi - |gg + m) 
~ Z^ Z^ ^(-«;0,0; gi,g 2 ) ^ 

91=0 <?2=0 

(7.25) 

fe =o i 1= o h=o V- + ^~ K > + s qi + 3^ - m) k 

J j 1 J rqi J r l 1=ni \J2- J j 2 +q 2 +h=n2 

Keeping the first K + 1 singular terms on the right and solving for the .^-coefficient with 
qi = ni, = Ti2, we have 
E(—k; 0,0; ni, n 2 ) 



{ I 1 + + 3 n l + 3 n 2 — Wjfc 



X 

fc: 

— 7T m! 



sin(7r/i(ft)) r(l + + m)r(— k) — |ni — |n2 + m 

" !; r(-/i(-«)-| ?1 -i(ft + TO 



(2KSo) m 6(«; m, ni, 712) 



Y^ Y^ ztV n n \ V 1^ V / 3^1 3^2 t 

" ££ E( - KiMi "'* ) r(^(- (t )-|„ 1 -^ + m ) 

(gi,92)/(ni,n 2 ) 

''' ? (1 + /*("«) + §9i + §*>)/ 



1 • - ■ - A / 1 — ^ 



X 



x f — (-2Kt 10 ) J1 J I — (-2^ 20 ) J2 , 

/ j 1 +q 1 +l 1 =n 1 \^2- / j 2 +q2+h=ri2 

(7.26) 

This is essentially (4.15) because of (4.11) and completes the proof of Theorem 1. 

7. 5 Choice of the computational parameter A 

If D s = D 6 = 0, the recurrence relation (2.7) for b(n; m, 0, 0) reduces to a two-term 
relation, and we obtain 

6(k; m, 0, 0) = (2KS )~ m — 3 V , 3 3 V — . (7.27) 

ml 
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Rewriting (4.16) by means of the identity 

(x) k -j __ (x) k (-k). 



(k-j)\ k\ (l-x-k) 3 

we then have 



(7.28) 



H( ., M , 0,0,0)- fc! ^ (l-^ ) - % (l + M -«)), j ! • ^ 

Because of 3^(k) = A + t(k) — 3 according to (2.5) and (2.11) and t(k) = 3 — t(—k) 
according to (1.6), the series is a terminating one-balanced hypergeometric series at unit 
argument, which can be summed by the theorem of Saalschiitz [11], so that 

H(- K ;t,0,0;0,0)= ( "» £ . 1 + J A > t( ^+,^ )t . (7.30) 

(1 + fj,{-K)) k k\ 

In total, we get 

-7T r(-iL+ir(K)+m)r(|L+ir(K) + m) 



e(-K;0,0) = 



r(-|L + ir(K))r(iL + \t(k)) TUX + ±t(k) + m)T(-\X + ±t(k) + m) 



3 ' 

>< 1 + e ( ;A L j #4 L + It? + q^- 1 ) 1 

(|A+|r(-«)-m) fc A;! 

(7.31) 

For A = — L or A = L, the terms with k = 1, 2, . . . of the asymptotic series all vanish and 
the factor in front of the series becomes independent of m. More generally this means 
that some terms which might become quite large when L is not small can be removed by 
such a choice of A from the asymptotic series and incorporated in the m-dependence of the 
function in front of the series. 



31 



References 

[I] W. Balser, W. B. Jurkat, and D. A. Lutz, Transfer of connection problems for mero- 
morphic differential equations of rank > 2 and representations of solutions, J. Math. Anal. 
Appl. 85 (1982) 488-542. 

[2] B. L. J. Braaksma, Multisummability and Stokes multipliers of linear meromorphic 
differential equations, J. Differ. Equations 92 (1991) 45-75. 

[3] W. Biihring, The characteristic exponent of second-order linear differential equations 
with two irregular singular points, Proc. Amer. Math. Soc. 118 (1993) 801-812. 

[4] A. Duval, Triconfluent Heun's equation, in: A. Ronveaux, ed., Heun's differential 
equations (Oxford University Press, Oxford, New York, Tokyo, 1995). 

[5] R. Hoeppner and R. Schafke, On the remainders of asymptotic expansions of solutions 
of linear differential equations near irregular singular points of higher rank, Math. Nachr. 
205 (1999) 89-113. 

[6] G. K. Immink, A note on the relationship between Stokes multipliers and formal 
solutions of analytic differential equations, SIAM J. Math. Anal. 21 (1990) 782-792. 

[7] E. L. Ince, Ordinary differential equations ( Dover, New York, 1956). 

[8] M. Kohno, Stokes phenomenon for general linear ordinary differential equations with 
two singular points, in: Japan - United States seminar on ordinary differential and func- 
tional equations, Lecture Notes in Mathematics 243 (Springer, 1971) 310-314. 

[9] M. Kohno, A two point connection problem for general linear ordinary differential 
equations, Hiroshima Math. J. 4 (1974) 293-338. 

[10] M. Loday-Richaud, Calcul des invariants de Birkhoff des systemes d'ordre deux, 
Funk. Ekvac. 33 (1990) 161-225. 

[II] Y. L. Luke, The special functions and their approximations, Vol. 1 (Academic Press, 
New York, 1969). 

[12] J. Martinet and J.-P. Ramis, Elementary acceleration and multisummability I, Ann. 
Inst. Henri Poincare, Phys. Theor. 54 (1991) 331-401. 

[13] F. Naundorf, A connection problem for second order linear differential equations 
with two irregular singular points, SIAM J. Math. Anal. 7 (1976) 157-175. 

[14] F. Naundorf, Ein Verfahren zur Berechnung der charakteristischen Exponenten 

32 



von linearen Differentialgleichungen zweiter Ordnung mit zwei stark singularen Stellen, 
Z. Angew. Math. Mech. 57 (1977) 47-49. 

[15] A. B. Olde Daalhuis and F. W. J. Olver, On the calculation of Stokes multipliers for 
linear differential equations of the second order, Methods Appl. Anal. 2 (1995) 348-367. 

[16] F. W. J. Olver, Asymptotics and special functions (Academic Press, New York 1974). 

[17] R. Schafke and D. Schmidt, The connection problem for general linear ordinary 
differential equations at two regular singular points with applications in the theory of 
special functions, SI AM J. Math. Anal. 11 (1980) 848-862. 

[18] Y. Sibuya, Global theory of a second-order linear ordinary differential equation with 
a polynomial coefficient (North-Holland, Amsterdam, New York, 1975). 

[19] V. S. Varadarajan, Linear meromorphic differential equations: a modern point of 
view, Bull. Amer. Math. Soc. 33 (1996) 1-42. 

[20] E. Wagenfiihrer, Die Determinantenmethode zur Berechnung des charakteristischen 
Exponenten der endlichen Hillschen Differentialgleichung, Numer. Math. 35 (1980) 405- 
420. 

[21] E. Wagenfiihrer and H. Lang, Berechnung des charakteristischen Exponenten der 
endlichen Hillschen Differentialgleichung durch numerische Integration, Numer. Math. 32 
(1979) 31-50. 

[22] E. T. Whittaker and G. N. Watson, A course of modern analysis (Cambridge Uni- 
versity Press, London , 1927). 

Wolfgang Biihring 
Physikalisches Institut 
Universitat Heidelberg 
Philosophenweg 12 
69120 Heidelberg 
GERMANY 



buehring@physi.uni-heidelberg.de 



33 



